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ABSTRACT 

This paper investigates the efficiency of an alternative to ratio estimator 
under the super population model with uncorrelated errors and a gamma- 
distributed auxiliary variable. Comparisons with usual ratio and unbiased 
estimators are also made. 
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I. INTRODUCTION 

It is well known that the ratio method of estimation occupies an 
important place in sample surveys. When the study variate y and the 
auxiliary variate x is positively (high) correlated, the ratio method of 
estimation is quite effective in estimating the population mean of the 
study variate y utilizing the information on auxiliary variate x. 

Consider a finite population with N units and let Xi and yi denote 
the values for two positively correlated variates x and y respectively for 
the ith unit in this population, i=l,2,...,N. Assume that the population 
mean X of x is known. Eet x and y be the sample means of x and y 
respectively based on a simple random sample of size n (n < N) units 
drawn without replacement scheme. Then the classical ratio estimator for 
7 is defined by 



A=T(A/x) (1.1) 

The bias and mean square error (MSE) of y,are, up to second order 
moments, 

B{y^) = x(RS\-S^,)lX 

M{y^)=x[sl+R^Sl-2RS^X 
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( 1 . 2 ) 

(1.3) 




where X = {n - n)l{nN) , 

K=YlX, sl={N-iyf^{y,-Y)\sl=(^-iy X (x.-X)% 

;=1 /=1 



and5,,= (N-l)-‘£ (y.-7)(x,-X). 



It is clear from (1.3) that M(yJ will be minimum when 

R=s^Js^, =p , (1.4) 

where p is the regression coefficient of y on x. Also for R = /? , 
the bias of y, in ( 1.2) is zero. That is, is almost unbiased for Y . 

Let E ( y|x ) = a + p X be the line of regression of y on x , 

where E denotes averaging over all possible sample design simple 
random sampling without replacement (SRSWOR).Then p^S^JSl 

and Y = a + p X so that, in general , 

R = ( alX )+ p (1.5) 

It is obvious from (1.4) and (1.5) that any transformation that brings the 
ratio of population means closer to p will be helpful in reducing the 
mean square error (MSE) as well as the bias of the ratio estimator y^ . 
This led Srivenkataramana and Tracy (1986) to suggest an alternative to 
ratio estimator y^ as 

y„ =z(x/x)+T = y^ -T{( a/x)-1} (1-6) 

which is based on the transformation 

z=y-A , (1.7) 

where E(z) = Z(= 7 - 4) and A is a suitably chosen scalar. 



In this paper exact expressions of bias and MSE of y„ are worked 
out under a super population model and compared with the usual ratio 
estimator. 



2. THE SUPER POPULATION MODEL 
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Following Durbin (1959) and Rao (1968) it is assumed that the 
finite population under consideration is itself a random sample from a 
super population and the relation between x and y is of the form: 

y. =a + /?Xi+ Ui ; (i= 1,2,...,N) 

where a and /? are unknown real constants; m,’s are uncorrelated 
random errors with conditional (given Xi) expectations 

E {u^ \x - ) = 0 
E {u^ |x,. ) = (5 xf 

(i-l,2,....,N), o{5{co , o<g<2 and Xi are independently identically 
distributed ( i.i.d.) with a common gamma density 

G(^) = e-V'Vr^,x)o, 2<^(oo . (2.1) 

We will write Ex to denote expectation operator with respect to the 
common distribution of Xi (i=l,2,3,...,N) and Ex Eg, as the over all 
expectation operator for the model. We denote a design by p and the 
design expectation Ep, for instance, see Chaudhuri and Adhikary 
(1983,89) and Shah and Gupta (1987). Eet ‘s’ denote a simple random 
sample of N distinict labels chosen without replacement out of 
i= 1,2,3 N. Then 

X(=N A) = ^ Xi + ^ X, 

iGs i^s 

Following Rao and Webster (1966) we will utilize the 
distributional properties of Xj / Xi , ^x, , , ^x. / ^x,. in our 

IGS i^s IGS i^s 

subsequent derivations. 

3. THE BIAS AND MEAN SQUARE ERROR 
The estimator in (1.6) can be written as 
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based on a simple random sample of n distinct labels chosen without 
replacement out of i = 1,2,... ,N. 

The bias 

B-Ep(T,-7) (3.2) 

of Ta has model expectation Em(B) which works out as follows: 



Era ( B ( J J ) - Ep Ex Ec 
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ExEc ( a + j3X) 



= EpEx 
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-A {(«/A)(l + (A-«)^/(«0-l))-l}-a 
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(N-n) {a - A)I N{ne - l) 



-A [-{N-n)lN + {{N-n)neiN{ne-\)]] 

(3.3) 



For SRSWOR sampling scheme , the mean square error 

Mfe)=Ep(jj,-Ff (3.4) 

of Fa has the following formula for model expectations 
Em(M(Fj) : 

Em (m(f^ )) = [e^ (m(f^ ))+{N -n\Nn0 + 2N -InfA^ - 2 A a)l N\n0 - \\n 0 - 2)J 

(3.5) 

where 

M(F)=£>F-rf (3.6) 

is the MSB of under SRSWOR scheme has the model expectation 

EAM{y^)={{N-n)lN^] 

A{n0 + g- 1)(« 0 + g-2) + n0{N0 -n0 + \f\ 

\ {Nn0 + 2N -2n)a^ \ J J r(^ + g) 

I («^-l)(«^-2) I {n0 + g -\^n0 + g -2) T0 

(3.7) 

[See,Rao{\96%,pA29)] 

Further, we note that for SRSWOR sampling scheme, the bias 

5(a) = ^,(a-e) (3.8) 

of usual ratio estimator has the model expectation 

^m{B{y^)) = {N-n)al [n0-\) (3.9) 

We note from (3.3) and (3.9) that 

if 
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(3.10) 




which is the same as in (3.10). 

Thus we state the following theorem: 

Theorem 3.1 : The estimator is less biased as well as more efficient 
than usual ratio estimator y^ if 

o{A{2a {a ^ o) 



i . e . when A lies between o and 2a . 

Therefore , when intercept term a{^ o) in the model (2.1) is sizable , there 
will be sufficient flexibility in picking A. 

It is to be noted that for o,y^ is unbiased and efficient than y^ . 
The minimization of (3.5) with respect to A leads to 

A ^ a = Aopt(say) (3.12) 

Substitution of (3.12) in (3.5) yields the minimum value of 

EA^{yo))as 

min E (mIv {N -\)s[{ne + g-\\n6 + g-2)+ne{N6 -nO + \)]T{e + g) 

■ („^ + g_l)(„^ + g_2) 

(3.13) 

which equals to E^{M{y^)) when a = o. 

It is interesting to note that when A =a,y^ is unbiased and attained its 
minimum average MSE in model (2.1). 
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In practice the value of a will have to be assessed, at the estimation 
stage, to be used as A. To assess a, we may use scatter diagram of y 
versus x for data from a pilot study, or a part of the data from the actual 
study and judge the y-intercept of the best fitting line. 

From (3.7) and (3.13) we have 

(T. )) - min -E^ [M (y J) = {(iV - «)(iV« 0 + 2iV - 2« )a ^ }/{iV ^ (« 0 - 1)(« 0 - 2)} 

) o (3.14) 

which shows that is more efficient than ratio estimator when A=a 
is known exactly. For a = o 

min.E„, {M{yJ] = E^{M{y^)) 

For SRSWOR , the variance 

N{y) = EXy-fJ 

of usual unbiased estimator has the model expectation: 

[v{y)) = {N- n)[f}^0 + 1^(0 + g)lY0]]l nN (3.1 7) 

The expressions of E^{M[y^)) and E^{v{y)) are not easy task to compare 
algebraically. Therefore in order to facilitate the comparison, denoting 

-100^Jf(t ))/fJm(tJ) and F, = 100FjF(yJ)/FjM(yJ), 

we present below in tables 1,2,3, the values of the relative efficiencies of 

T^with respect to y and y^ for a few combination of the parametric 
values under the model (2.1). Values are given for N = 60 , 
0 = 2.0,0 = 8,a = 0.5, 1.0, 1.5, = 0.5, 1.0, 1.5 and g = 0.0, 0.5,1.0,1.5,2.0. 

The ranges of A, for y^ to be better than y^ for given a = 0.5, 1.0, 1.5 are 
respectively ( 0,1), ( 0,2), (0,3). This clearly indicates that as the size of 
a increases the range of A for y^ to be better than y^ increases i.e. 
flexibility of choosing A increases. 

We have made the following observations from the tables 1,2 and 3 : 

(i) As g increases both Ei and E 2 decrease. When n increases 
El increases while E 2 decreases. 



(3.15) 

(3.16) 
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(ii) As a increases ( i.e. if the intercept term a departs from 
origin in positive direction) relative efficiency of with 

respect to y decreases while E 2 increases. 

(iii) As p increases Ei increases for fixed g while E 2 is unaffected. 

(iv) The maximum gain in efficiency is observed over y as well as 
over A if A coincide with the value of a. Finally, the 
estimator y^is to be preferred when the intercept term a 
departs substantially from origin. 
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